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FIGURE CAPTIONS 


Figure 1. 

Figure 2. 

Figure 3. 
Figure 4. 
Figure 5. 
Figure 6. 


Isotherms for the Benard convection cell with 
P = 10. 

Isotherms for the Benard convection cell with 
P = 80, 

Driven cavity. 

Stream function contours at R = 100 (41 x 41 grid) . 
Stream function contours at R = 1000 (61 x 61 grid). 
Single vortex convected over a porous wall. 
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SECTION 1 


INTRODUCTION 


The turbulent bursting phenomenon is recognized as a 
dominant feature of time-dependent turbulent boundary-layer 
flow. The nature of the boundary layer formed on a flat plate 
located beneath a convecting rectilinear vortex embedded in a 
uniform flow was investigated by Doligalski and Walker (1) . 

They obtained numerical solutions for the temporal development 
of the boundary layer induced by the motion of a rectilinear 
vortex. The boundary layer was expected to erupt from the wall 
into the inviscid flow. 

In the present study the flat plate is replaced by a porous 
wall where the velocity at the wall obeys Darcy's law. The 
point of inflexion of the boundary layer velocity profile is 
controlled by varying the strength and/or the location of the 
vortex. Our interest is focused on how the inflected velocity 
profile responds to wall mass transfer induced by the motion 
of the vortex. 

The governing equations for a two-dimensional incompressible 
viscous flow are the vorticity transport equation and the stream 
function equation. These equations can be solved simultaneously 
at each time step by means of a finite difference scheme. We 
propose to use the weighted-mean scheme of Fiadeiro and Veronis 
(2) for the vorticity transport equation. This scheme is first 
applied to a series of test problems to determine its accuracy, 
stability and efficiency. For the stream function equation, we 
use the usual central difference scheme. 
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SECTION 2 - THE l-JEIGHTED-MEAN SCHEME 


The advection-dif fusion equation is written in conservative 

form : 
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where v is the diffusivity constant 

u and V are the velocity components in the 
X and y directions, respectively and 
oj is the vorticity 

For cases when the vorticity field is subject to strong advection, 
centered-f inite difference procedures are inefficient; consequently, 
the weighted-mean scheme is adopted for the present investigation. 

An explicit (leap-frog) centered-time difference in equation 
(1) is used except for the central term that results from 

the spatial differencing evaluated as 

where the superscript denotes the time step. When the weighted- 
mean scheme is applied to the spatial derivatives the equation 
for the calculation of given values at and is 
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where At is the time increment 

h is the spatial increment 
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with all a's evaluated at time step n, we have: 
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The error of approximation is 0(At , h ). The formulation 
uses only a five point operator for two-dimensional flow and 
is antisymmetric in relation to the velocity field. When the 
component change sign, the coefficients upstream and downstream 
of the point are autom.atically reversed, a feature particularly 
useful in computer programming because the sign of the velocity 
components is not generally known a priori. 
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SECTION 3 


TEST PROBLEM 1 


The one-dimensional advection-dif fusion problem presented 
by Roache (3) to illustrate the weakness of the standard central 
difference relation is as follows: 


i ^ + H 

V 3t V 9x 


9x^ 


= 0 


(3) 


v;ith boundary conditions 
o)=0 atx=0 
o)=l atx=l 

Regardless of the initial condition, the steady state 
solution for incompressible flow is 
ux/v 


_ 1 - e 
^ 


1 - e 

The time derivative is written as first-order forward-dif- 
ference approximation and the weighted-mean scheme is applied 
to the spatial derivatives. 

Steady state results appear in Table 1 for weighted-mean 
scheme, upwind differencing scheme and central difference 
scheme for various values of cell Reynolds number 
uh 


Rc = 


V 


The steady state form of equation (3) is 
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It can be shovm that the solution of the finite difference 
analog of this equation agrees identically at the points x^ 
with the solution of the differential equation (4). We may 
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observe this remarkable result in Table 1. Table 1 also reveals 
that the central difference scheme produces wiggles for values 
of Rc greater than 2.0. The upwind differencing scheme, although 
stable, is considerably in error at all values of Rc. 
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SECTION 4 


TEST PROBLEM 2 


The problem is to determine the temperature field in a 
two-dimensional square, Benard convection cell for a fluid 
heated uniformly from below and cooled from above. The length 
scale L and the maximum amplitude U of velocity are used in 
the non-dimensionalization . 

The non-dimensional temperature equation has the form: 

„ ,9(uT) , 9(vT) , _ 9^T . 9^T 

where P = Peclet number = ^ 

K 

K = thermal conductivity 
The dimensionless velocity field is given by 

u = - sin X cos y v = cos x sin y 0 £x£tt 

0<y<Ti 

Boundary conditions on T are 
T = 1 at y = 0 
T = 0 at y = TT 

=0 at X = 0 ,TT 

dX 

A relatively coarse grid system of 25 x 25 mesh points was 
used for different values of P. Although the scheme assures 
convergence to all iterative methods, we would like to compare 
the rates of convergences in Table 2. Line SOR combined with 
multigrid techniques are approximately four times faster than 
ADI method. Figures 1 and 2 show isotherms of the fluid for 
Peclet numbers 10 and 80. The intertwining tongue for the 
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case P = 80 indicates the effect of strong convection. A final 
point that we would like to bring out for the strongly con- 
vective case (P = 80) is that distortions of isotherms were not 
introduced with relatively coarse grid (25 x 25) . 
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SECTION 5 


TEST PROBLEM 3 


We consider the flow in a driven square cavity. This problem 
is a standard test case for evaluating the accuracy, stability 
and efficiency of numerical schemes. For simplicity we choose 
a square with sides of length L = 1. The upper surface (fig. 

3) moves to the right with a constant transverse velocity U = 1. 

The flow of a viscous incompressible fluid in a square 
cavity is governed by the following coupled equations: 
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where ij; = stream function 
u) = vorticity 
R = Reynolds number - 
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The vorticity transport equation (5) becomes equation (2) in 
finite difference approximations. The stream function equation 
(6) can be discretized using central difference formula. For 
concreteness, we will describe the computation cycle for the 
time-dependent approach. 

The calculation starts at t = 0 where and u) are known 
everywhere. From equation (5) we obtain the vorticity u for 
all interiors points at t = At. The numerical solution of (6) 
gives at t = At and hence u and v. The last step of the 
computational cycle is to update the boundary values of co 
using the most current values of i> and co at the interior points. 
This computational cycle is repeated until the steady state 
is reached to a specified convergence level. If we use the 
steady state form of equation (5) , instead of "marching" the 
vorticity transport forward in time, new values of a; are obtained 
by an iterative process. The way in which the boundary values 
of the vorticity are approximated affects both the rate of 
convergence and the accuracy of the solutions. This fact has 
been noted by many authors. In the present case we use a second 
order form for evaluating wall vorticity given by 
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Studies indicated that instead of using formula (7) at time step 


ri , the average (ui 
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wall 


+ (jO 


n-1 

wall 


)/2 yielded a more stable process 


for the time-dependent and steady state approaches. Contour 
plots of stream function are shown in figure 4 and 5 for various 




5-2 


grid sizes and Reynolds number. Numerical values also are listed 
in Table 3 for comparison. 
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SECTION 6 


COMMENTS ON THE VORTEX PROBLEM 


The problem is shown schematically in figure 6. An isolated 
vortex embedded inside a boundary layer convects to the right 
over a porous wall. The velocity at the wall is given by 
Darcy's law: 


where c is a constant 

Pr is a reference pressure 

The wall mass transfer mechanism is driven by the wall 
pressure gradient induced by the motion of the vortex. Consequently, 
it is of interest to examine the effect of wall pressure 
fluctuations on velocity profile with an inflexion point inside 
the boundary layer. The point of inflexion initially may be 
created by varying the position and/or the strength of the 
vortex . 

It is worthwhile to first consider a non-porous wall. 

The effect of two counter-rotating vortices of opposite and 
equal strength is equivalent to having an impermeable wall 
between them. In an inviscid flow, a dimensionless stream 
function of the vortices may be defined by 


V (x,t) = c(Pr- P (x,t)) 

W W 1 



2 



2 




where F is the strength of the vortex 


r is the core radius 
c 


(x ,y ) is the center of the vortex 
o o 
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Let the vortex be convected in an inviscid flow. The convection 

speed of the vortex is composed of 2 parts consisting of the local 

mean velocity and the self induced velocity due to the vortex 

image below the wall. Note that the inviscid solution predicts 

that the vortex will remain at constant heiaht y . However, 

-^o ' 

the inviscid solution is not valid for our realistic problem 
because it fails to satisfy the no-slip boundary condition on 
the wall. 

We have also tested the possibility of satisfying the no- 
slip condition by means of a smoothing function, s. The modified 
stream function is defined as: 

= s ij; 

modif orig 

and s vanishes strongly near the wall, to ensure zero derivatives 

of i' , ■ j. (zero velocity). Far from the wall, s is identically 
modi t 

1. However, this new stream, function produced a reverse flow 
near the wall, introducing new complications in the flow pattern. 
It was therefore abandoned. 


SECTION 7 


CONCLUDING REMARKS 


The weighted-mean scheme has several significant advantages: 
it is locally a conservative scheme neither creating nor 
destroying the advected quantity artifically. ' 

it is stable for all grid spacings or cell Reynolds numbers. 

The scheme thus has the practical advantage of yielding 
solutions for relatively large values of grid size and is, 
therefore, economical with computer time. One can use this 
crude solution to construct a good initial guess solution 
for a finer grid. 

it becomes a central difference scheme for strongly dif- 
fusive cases and an upwind difference scheme for strongly 
advective cases. Furthermore, when the components of the 
velocity change sign, the coefficients upstream and downstream 
of the point are automatically reversed, a feature particularly 
useful in computer programming because the sign of the 
velocity components need not be known a priori. 

The only major disadvantage is the loss in speed due to 
the computation of the hyperbolic tangent at each grid point 
in m.ultidimensional problems. 

For the vortex problem, the initial conditions need to 
correspond to some real initial situation since we are interested 
in the transient solution. Therefore one needs to satisfy 
initially the no-slip condition. One suggestion is that the 
vortex at time t equals zero, should be far away from the wall 
in the inviscid region so that the velocity induced by the 
vortex vanishes at the wall. 
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Table 1 


Rc = 1.0 
Grid point 



1 

2 

3 

4 

5 

6 

ES 

0.0 

0.0117 

0.0433 

0.1295 

0.3636 

1.0 

WMS 

o 

o 

0.0117 

0.0433 

0.1295 

0.3636 

1 . 0 

UDS 

O 

o 

0.0323 

0.0968 

0.2258 

0.4839 

1.0 

CDS 

0 . 0 

0.0083 

0.0331 

Rc = 2.0 
Grid point 

0.1074 

0.3306 

1 . 0 


1 

2 

3 

4 

5 

6 

ES 

O 

o 

0.0003 

0.0024 

0.0183 

0.1353 

1 . 0 

WMS 

o 

o 

0.0003 

0.0024 

0.0183 

0.1353 

1.0 

UDS 

0 . 0 

0 . 0082 

0.0331 

0.1074 

0.3306 

1.0 

CDS 

o 

o 

0.0000 

0.0000 

Rc = 4.0 
Grid point 

0.0000 

0.0000 

1 . 0 


1 

2 

3 

4 

5 

6 

ES 

o 

o 

0.0000 

0 . 0000 

0 . 0003 

0.0183 

1.0 

WMS 

o 

o 

0.0000 

0.0000 

0.0003 

0.0183 

O 
1 — 1 

UDS 

o 

o 

0 . 0013 

0.0077 

0.0397 

0.1997 

o 
1 — 1 

CDS 

0 . 0 

0.0164 

-.0328 

0.1148 

-0 . 3279 

o 
1 — 1 


ES = exact solution 
WMS = weighted-mean scheme 

= upwind differencing scheme 
= central difference scheme 


UDS 

CDS 


Table 2 


Benard cell (25 x 25) 


-MULTI GRID (8) 


p 

SOR 

ADI 

in work unit (wu) 




SOR 

Line SOR 

10 

140 

85 

38 

20 

20 

188 

110 

35 

27 

40 

A 

140 

44 

38 

80 

A 

205 

83 

51 

160 

A 

A 

u- 

101 

96 

A The 
not 

results for those 
coinputed . 

specific 

Peclet numbers 

were 


Table 3 


R 

Reference 

Grid 

VC 

Iterations 

100 

h’MS 

21 X 

21 

-0.094 

55 (TD) 


WMS 

41 

41 

-0.103 

120 (ss) 


Burggraf (4) 

41 

41 

-0.101 



B + D (5) 

51 

51 

-0.103 



N + KP (6) 

51 

51 

-0.102 


1000 

WMS 

61 

61 

-0.092 

210 (ss) 


B + D (5) 

51 

51 

-0.081 



N + KP (6) 

51 

51 

-0.097 



Gosman et al. (7) 

81 

81 

-0.099 



WMS = V7eighted-mean scheme 

TD = time dependent approach 

ss = steady state approach 

t = stream function at vortex center 
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